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Abstract—Existence of symmetric waves in open metal- 
dielectric waveguides, a dielectric rod and the Goubau 
line, is proven by analyzing the of functional proper- 
ties of the dispersion equations (DEs) and parameter- 
differentiation method, applied to the analytical and 
numerical solution of the DEs. Various limiting cases 
are investigated. Reduction to singular Sturm-—Liouville 
boundary eigenvalue problems on the half-line is per- 
formed. Principal and higher-order surface waves are 
investigated. 


Index Terms—Goubau line, Surface wave, Disper- 


sion equation. 


I. INTRODUCTION 


Analysis of wave propagation in open metal—dielectric 
waveguides constitutes an important class of vector elec- 
tromagnetic problems. A conducting cylinder covered by a 
concentric dielectric layer (see Fig. 1) is the simplest type 
of such guiding structures. The first studies date back to A. 
Sommerfeld who wrote that, ” Taking into account the finite 
conductivity of a straight, round wire, a possible solution of 
Maxwell’s equations is a surface wave that travels along the 
wire and is non-radiating. This will be a TM wave at least 
consisting of the fundamental mode” [1]. Harms [2] was the 
first to obtain a dispersion equation (DE) for a conducting 
wire surrounded by a dielectric sheath. This structure was 
later called the Goubau line (GL) [3]. A thin dielectric cover 
in the Goubau line plays the role of conductivity and this 
line allows to transmit electromagnetic waves along a single 
wire when non-radiating surface waves are bound to the wire. 
Since the early 1950s, a wide variety of studies and applica- 
tions of Goubau lines have been performed and published. 
We note, for instance, paper [4] dealing with theoretical and 
experimental investigations of surface waves along coated 
and corrugated conductors and guides with capacitive surface 
impedance, and book [5] which summarizes some aspects of 
the theory of surface waves. A review of different types of 
Goubau-type surface-wave transmission lines and a summary 
of the earlier results of Sommerfeld and Harms-Goubau lines 
are presented in [6], as well as the DEs for surface waves 
and various known applications. Mathematical and numerical 
investigations of the mode spectrum of Goubau lines are 
performed in [7-9]. A comprehensive list of references of 
more recent results can be found e.g., in [10, 11]. Several 
recent applications of Goubau lines are reported in [12, 13]. 
In spite of a long history of the studies partially reviewed 


above, many basic problems that have been solved for empty 
shielded waveguides (existence of normal waves, discrete- 
ness and localization of their spectrum on the complex plane, 
and so on) still remain unsolved for open metallo—dielectric 
waveguides and even for the Goubau line. In particular, to the 
best of our knowledge a satisfactory mathematical analysis of 
multi-parameter DEs has yet to be carried out. This fact has 
become a driving force for us to carry out a systematic study 
of the corresponding mathematical issues using the recently 
obtained fundamental results in the mathematical theory of 
electromagnetic wave propagation [14, 15]. In these works, 
it is shown in particular that the spectrum of a broad family 
of metal-dielectric waveguides consists of isolated points on 
the complex plane (continuous spectrum is absent) and is 
localized in a strip; the analysis is based on the operator 
spectral theory methods. 


In this paper we consider two basic structures: a dielectric 
waveguide of circular cross-section and the Goubau line. 
We reduce the problems of normal waves to the singular 
Sturm—Liouville boundary eigenvalue problems on the half- 
line with a discontinuous (piecewise constant) coefficient. It 
should be noted that the boundary conditions in the latter 
problem are complex functions of the spectral parameter 
and therefore there are no results concerning the spectrum 
existence and distribution in the classical literature on Sturm— 
Liouville boundary value problems. In such cases, when the 
boundary operator depends in a complicated manner on the 
complex spectral parameter, each Sturm—Liouville boundary 
eigenvalue value problem or a family of such problems 
should be investigated separately using the analysis of this 
particular dependence. The determination of eigenvalues is 
reduced to certain nonlinear equations F'(z,@) = 0, called 
the DE for normal waves. Here, z is the spectral parameter 
(the longitudinal wavenumber of the running wave) and 
æ = (a1, 42,...) is the vector of nonspectral parameters of 
a particular problem. DE specifies z = z(@) as an implicit 
function. In many practical examples, only one nonspectral 
parameter is varied. We prove the existence of roots of DE 
by reducing it to the form F(z, a1) = F2(z), where function 
F(z) is characterized by infinite number of singular points 
and takes all real values on the intervals between these points, 
and the left-hand side is a continuous function of a fixed 
sign on a certain interval. Using the method of the parameter 
differentiation, DEs are reduced to the Cauchy problems to 
find an implicit function z = z(a,). The determination of 
the propagation constants of symmetric waves propagating 
in the Goubau line with a thin, perfectly conducting cylinder 
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placed inside a dielectric rod and with a thin dielectric layer 
are considered when the radius of the cylinder (or the width 
of the thin layer) is taken as a small parameter (a,). In 
its mathematical completeness including the development 
of the Sturm-Liouville theory, the present study is the first 
and necessary step towards the analysis of dielectric and 
metal—dielectric waveguides (fibers, Goubau lines) including 
complex waves, resonance and interaction phenomena, and 
structures with more complicated coating formed by linear 
and nonlinear layered or radially inhomogeneous media. 


$ 


p 


Fig. 1: The Goubau line, where a and b are the radii of the internal (perfectly 
conducting) and external (dielectric) cylinders. 





II. STATEMENT 


Consider the propagation of symmetric eigenwaves, de- 
scribed in terms of nontrivial solutions to homogeneous 
Maxwell’s equations, in the dielectric rod with a radius a and 
the Goubau line, in which a and b, b > a, are the radii of 
the internal (perfectly conducting) and external (dielectric) 
cylinders. Longitudial wavenumbers k, of the symmetric 
waves are determined as 


k =p, Sa 
2 __ O ’ 5 
ko = ck? =f, Hea (dielectic rod), 
k — B87, r>b 
2 — 0 ’ ’ 
k? = l ch? —B2, a<r<b (Goubau line). 


Here, 6 is the wave propagation constant (spectral param- 
eter), € > 1 is relative permittivity of the homogeneous 
dielectric, and ko is the free-space wavenumber. 

We consider symmetric azimuthally-independent waves, 
having the nonzero components 


H = (0, Ho(r,z),0], E = [Fi(r,z),0, E3(r,z)], 0) 


ibd l 
E = 18 pe E3 = b(r)e—*, 
kg dr 
(2) 
iwe dò _ ig, 
H> = —— —e 5; 
k2 dr 


where @ is a cylindrical function satisfying the Bessel equa- 


tion: 
_ ld, do 24 


Determination of surface waves can be reduced to the fol- 
lowing singular Sturm—Liouville boundary eigenvalue prob- 
lems on the half-line with a discontinuous (piecewise con- 
stant) coefficient in the differential equation: 








_ _ | Ed z 
Lo = 0, r > 0, [Pll =a = $ J ni (4) 
d(r) +0, r > œ, 
$ € C![0, +20) N C? (0, a) N C? (a, +00) (5) 
for the dielectric rod; and 
Lo = 0, r>a, o(a) — 0, lln = 
e dd (6) 
= Fed BS = 0,¢(r) => 0, r > œ, 
b € C'[a, +00) N C° (a, b) N C? (b, +00) (7) 


for the Goubau line. 

The conditions at infinity in (4) and (6) and the form (2) 
can be taken as a definition of symmetric (surface) waves 
in open waveguides under study. The surface waves are 
described in terms of real-valued quantities; in particular, 
the boundary operators in (4) and (6) are defined (as real- 
valued functions of a real variable y or \ = 7°) on a certain 
interval 7. The spectrum of surface waves may be empty, 
or they may consist of several (real) points located on this 
interval. The results of the classical Sturm—Liouville theory 
concerning existence and distribution of the (real) spectrum 
are not applicable here, and we will perform in this work the 
corresponding analysis from the very beginning by reducing 
the boundary eigenvalue problems under study to DEs. This 
will, in turn, enable us to prove all the required statements 
and determine the eigenvalues numerically or analytically. 

Make the following designations: 


b 
=É, x = koay € — 7°, u = koave — 1, s=- >1. 
i a 

(3) 


For surface waves that must decay at infinity, potential 
function $(r) is represented in terms of cylindrical functions 
that tend to zero at infinity, so that 


a PP 
i= oa Th SG. (9) 
BJA T) Taa 
for the dielectric rod, and 
r >b, 


o(r) = 





AKo(LVu? — 2”), 
AG | Jo(Lx)¥o(x) — Yo(Zx)Jo(x)|,a < r < b 
(10) 
for the Goubau line. 
Here, J, (x), Yp(x), Ko(x), and K(x) denote, respec- 
tively, the Bessel and Neumann functions of the order 
k = 0,1 and the McDonald functions. 
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IHI. DIELECTRIC ROD 


By applying to (9) the condition of continuity of @ and its 
derivative at the point where the coefficient is discontinuous, 
we obtain the desired DE in the form [18] 





Jal) ___ z; Ho iV? = x?) 
*Ti(2) a H (ivu — 72) en 


which will be considered in what follows with respect to 
new variable x (instead of 3). Here, HP (.) denotes the 
zero-order Hankel function of the first kind. 

We rewrite (11) as 


Galz, u) = Hale), (12) 
Fp(#,u)=0, Fp(#,u) = Ga(a,u)— Fa(x), (13) 
where 
Ko(Vu2 — x?) 
= ~e\/u? — 2 14 
Ga(x, u) EVU — T TE (14) 
Taa) = zeot) coe) = re (15) 


In view of numerical solution and analysis of the limiting 
cases a — 0 and e — 1, it is reasonable to consider equation 
(12) with respect to the new independent variable 


2_ 72 = \/7?-1>0 
with 


E (0, Ve — 1) 


(Ca) =U 


y= koe 


where e€ > 1. 


40 








Fig. 2: Graphs of Gg(x, u) (——) and Fy(ax) (- - -) against x for € = 5, 
koa = 2.8, and u = 5.6 illustrating (graphical) solution of DE (12) for the 
dielectric rod. 


One may call function cotj(x) the cylindrical cotangent, 
because it exhibits all basic properties of the trigonometric 
cotangent; namely, cotj(x) takes all real values on every 
interval (uf, re), where ust) are zeros of cylindri- 


cal functions J(x) (k = 0,1), and cotj (x) > Ow € 


k)  (k+1 
(WO) wETY) b= 1,2)... 


Equation (13) determines x = x(u) as an implicit function 
of the parameter u > 0, and can be used to prove that this 
function exists for all u > u*, where u* > a Namely, 
using the properties of the cylindrical cotangent and function 
Ga(x,u), we arrive at the known result [2], which is: 


If 
u = koav e — 1 > po), 
then there exists a root 


* 0 1 
z= 2" € (uy, Hy’) 


of DE (13). If 
0 1 
ea a; 


then (13) has only one root. 


IV. GOUBAU LINE 
A. Dispersion equation 


By applying the transmission (continuity) condition to (10) 
we obtain the desired DE for surface waves in the Goubau 
line in terms of variables (8), which reads 


Tomus) zG ana] has) =: (16) 
where 
7 z a Kols Vu = a) 
G(x, U, 5) = Ey u? — z a (17) 
Fea 8 | pes =, (18) 
®1(x, 8) = Jo(sx)Yo(x) — Jo(x)Yo(sx), 
O5(x, s) = Jo(x)Y1 (sx) — Jı (sx)Yo(z). 


B. Principal surface wave 


Let us first prove the existence of the principal (funda- 
mental) surface wave for all u > 0 and s > 0. We denote by 

ye ) the zeros of the Neumann functions Yp(x) (k = 0,1). 
ae to [17], the zeros of the Bessel and Neumann 
functions alternate: 


(1) (0) 


1 
ps 02 024, ii 


1 0 
Sy ei a 
(1) (0) 


Sy ay, < Ha <. 


Denote by h? (q), and ht? (q) (m = 1,2,.. 
products of cylindrical functions 


Zoo(q, y) = Joly) Yo(ay) — Yo(y)Jo(ay), 
Zıolq, y) = Jı (qy)Yoly) — Yı (qy) Jo(y), 
For a fixed q > 1, all these zeros are real and simple points 
[17]; the sequences {h? (q)} and {h}0(q)} are increasing 
w.r.t the index m, and the zeros of their terms alternate on 
the positive semi-axis. The h!9(q) are decreasing functions 
of q, for q > 1. 
Thus, F,(x,s) in (18), can be represented as 


Zoo(s, £) 
Z10(S, x) 


(19) 


.) zeros of the 


(20) 


F,(a,8) = xcotj (x), cotj (x) = (21) 
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Ggix,u, 1.78) 
Ggix,u, 3.4) 10+ 
Gefx,u,4) 








0 1 2 3 4 
x 
a) Graphs of Gg(x, u, s) show that Gg is a slowly varying 
quantity w.r.t s and u and that Gg(x = u,u,s) = 0 Ys > 1. 
2 | ; 





b) Graphs of F(x, s) show that increasing s shifts the curves 


of cylindrical functions involved in the DE to the left. 


Fig. 3: Behavior of functions Ga(x, u, s) (17) and Fa(x, s) (18) for the 
Goubau line (e = 5) for different value s = 1.78 (— —), s = 3.4 (.....), and 
s = 4 (—). 


meme 





Fig. 4: Graphs of Zoo(q,y) and Z10(q, y) at € = 5 and q = 3.4 


where cot) (x) also resembles the cotangent and cylindrical 


cotangent. 
The domain of function F(x, u, s) is 


Dpr, lu, s): {x£ :|x|<u\Haelu,s)}, w>0, 
Hg(u, s) = {ht (8) :0 < h?! (8) Se= eee eres We 


Sl. 


Hg(u, s) = 0 (N(u) = 0) 
if all A°l(s) >u (mM =1,2,...N(u)). 


According to the conditions on parameters we will consider 
Fo(a, u,s) on the set 


Dy, (us 8) tts). 0 << Ue O > 0.8 > Lh. 


Feg(a,u, s) is a continuously differentiable function of x on 
DE. (u, s) for all u > 0 and s > 1. Below we will prove that 
N (u) > 0 specifies the number of surface waves propagating 
in the Goubau line for the given u > 0. 

From the explicit form of functions in (16) it follows that 
for all u > 0 and s > 1: (i) F,(x,s) takes, as a function of 
x, all positive values on the interval I(s) = (0, ht°(s)); and 
(ii) G(x, u, s) takes, as a function of x, all values from the 
interval Ir(s, u) = (0,T(s,u)) where 


Ko(su) 


LS 0) = U 





Consequently, there exists (at least) one root x° = x°(s, u) € 
I(s) of equation (16). This root corresponds to the funda- 
mental surface wave, which exists for all u > 0. The analysis 
is Clearly illustrated in Figs. 5 and 6, where the x-coordinate 
of the point of intersection on the interval between the origin 
and the first vertical asymptote is the root x° corresponding 
to the fundamental (no-cutoff) surface wave. 

Check the location of roots for large values of s. Using 
explicit form (20), one can show that, for s > 1, there exist 
two zeros of ®ə(x, s) when 


ge i yo) r 


In particular, 


ys? yo? 
s > 4% = 4.333: N=2: s> 3 =7.89. 
(0) (0) 
Yi Yi 
However, using the properties of the function that enter (17) 
and ( 18), one cannot say anything about the existence of 
(0) 


roots in the interval (0, #—). I will consider this interval 


below by reducing the examination to the analysis of the 
Cauchy problem. 

Sufficient conditions for the existence of two zeros of 
(x, s) can be obtained in a similar manner 


N 





s < y. (23) 
Thus, using (19) and (23), we find that when 
(0) (1) 
— fa a 2 
Yi Yi (24) 


= max{ 7.930044, 6.0763} = 7.930044, 
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the zeros of the numerator and denominator in (18) are 
disposed as follows: 


0< gen 2 gm, gor < r3 


Applying the properties of function cot) ( (x) similar to those 
of cot x (zeros and singular points of cotj (x) alternate, each 
zero is followed by a singular point, and the function takes 
all real values between the neighboring singular points), one 
can show that for 
yy” (0) 
u > a or kobve-— 1 >y’ = 3.957678 


and under the condition (24); that is, when 
kob 
7.930044 ’ 


there exists at least one root x = x, of equation (16) located 
in the interval (xf°", x”). 


koa < 


C. Higher-order waves, limiting cases and the parameter 
differentiation method 

The analysis of Section 2.3 is not applicable when s —> 
1. Therefore, we present an alternative approach based on 
parameter differentiation and implicit-function methods that 
enables one to consider various limiting cases. 

In (16), we introduce the change of variables 


ST = Y, 


If b is fixed, then t = 
(16) in new variables: 


Fo(y,t) = Gg(y, v) 


— 0 if a > 0. Rewrite equation 


œj 


— F,(y, t) =0, (25) 
where 
Ko( yv? — y?) 


G,(y,v) = e/V — y? 
g (¥ ) ER TE 


~ P, y, t 
A E N 
(y, t) 


Pı (y, t) = Jo(y)Yo(ty) — Jo(ty)Yo(y), 
Pa(y, t) = Jo(ty)Yi (y) — Ji(y)Yo(ty). 


One can show that for t — 0, when the Goubau line turns 
into a dielectric rod, equation (16) in the form (25) turns into 
DE (12) for symmetric waves of the dielectric rod. 

The roots y = y(t) = kob\/e — y? of (25) can be found 
as solutions to the Cauchy problem 


dy == P ( N ƏFa(y,t) 
a(y,t), hni Paa 26 
te eo awa = 
y 





where yo € (0, koby€) is any solution of equation (12). 
Taking into account the recurrence relations for cylindrical 
functions and the asymptotic behavior Yo(z) ~ Żlnz as 
z — 0 one can show that 


OF c(y,t) 





Oy y=yo, t=0 = Qyo), 
aG a(yo,0) Jo(yo) 2 Sit (yo) 
Quo) = + Jı (yo) Yo TETE AT 


The Cauchy problem (26) is solvable under the condition 


Q (yo) # 0, 


which provides also that ®g(y,t) is continuous in the 
vicinity of the point (y,t) = (yo, 0). 

Thus placing an infinitely thin (a ~ 0), perfectly con- 
ducting rod into a dielectric fiber (a = 0) creates a regular 
perturbation for the spectrum of eigenwaves of the dielectric 
fiber (rod). 

Let us examine another limiting case b— a — +0, or s > 
+1, when the Goubau line turns into a perfectly conducting 
cylinder with radius b. In the limit s — +1, (16) takes the 
form 
Ko(v u? — x?) 
Ky (Vu2 — x?) 


HO? (y fo SUP) =O) 


Hence, the propagation constants 


(ys 
zm 
=: (5) ; (27) 


z z) = 


u2 — r? — 0. 


or 


e 





— 0 where z4 ) denotes a 
OD 


are complex. Here, H 
(complex) zero of the Hankel function H iS Mz ) with Szm 
OSS Nan 

Other solutions in this limiting case can be obtained from 
the equation 


r =u = y=1. 


Implicit function « = x(s) can be determined as a solution 
to the Cauchy problem 


OF G(a,s) 





dx 
= Va(z, s), ~ a 
ds — —— OS 
l z(1) = Sii Va(z, s) T ’ (28) 


which is solvable under the condition 
OF g(a, s) 


Ox Zo 





TU; s= 
and We(x,s) is continuous in the vicinity of the point 
(x, s) = (u,1). We have 
OF g(a, s) 
Os 
OF g(a, s) 
Ox 





= Qi; 
Qı ~ ze [2s Inte ue a o(1)| >-0, £>-U. 


So, the Cauchy problem (28) is not solvable in the vicinity 
of the point (x, s) = (u,1). 

Thus, again we see that, generally, introduction of a 
thin dielectric layer is a not regular perturbation for the 
spectrum of eigenwaves produced by TEM wave (y = 1) or 
a decaying wave that corresponds to complex propagation 
constant yt) Below we will overcome this difficulty by 
replacing a singular form of the DE given in (25), by a 
regular one. 
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In order to prove the existence of a solution £x = x? = 


x°(s) to DE (16) in the vicinity of the point x = u 
corresponding to the principal surface wave consider the 
limiting case s — +1 and rewrite equation (16) in the regular 
form 


Fo(a, u, s) =G (x,u, 8) ®o(a, s) — 


29 
— zı (z,s)=0, O<z<u. oe 


We have Fo(x = u, u, 1) = 0 for all u > 0 and OF olthsthst) A 
0 (e.g. EES ~ 3.6538 at u = 3.999... and s = 1). 
Consequently, x = x(s) = x? (s) is a regular perturbation of 
the value x = u (or y = 1) attained at s = 1 and can be 
determined as a solution to the Cauchy problem 


0 ~ OFo(a,u,s) 
S Pets), panei ee 
a(1) =u. ile: OFo(w,u,s) ` 

Ox 


Note that here, u is a fixed positive number. 

Let us now consider the higher-order surface waves. Note 
first of all that: (1) for every u > 0 there exists a (minimal 
positive) zero h? (s) of function ®,(z,s) such that 

u =h? (s?) for one and only one 


(30) 


sap Ss) Si: IWS 12a 


(ii) according to the explicit form (29) of the DE Folz = 
u, U, s? ) = 0 for any u > 0 and 2 Fouts) Æ 0. Numbers 


Sņ can be determined, for every given u > 0, by numerically 


solving the equation 
Pı (u, s) = 0. (31) 


Thus z = z™ (s) = x™ (s) can be determined as a regular 
perturbation of the value x = h?, attained at s = s9,: for a 
fixed u > 0, by finding the implicit function 7” = xz™ (s) 
as a (real) solution to the Cauchy problem similar to (28) 


drt” _ T 0 
Te LOS), BSS. 





GS So aa De oa 


Note that x” (s) satisfies the condition (see Fig. 8) 
h? (s) > £™ (s) > h?! (s), m=1,2,...No(u). (33) 
We can verify the solvability conditions: 


OF (2, u, s) 


Ox PY 


z=h9, (5), s=89, 





and F(a, u, s) is continuous in the vicinity of point A” = 
(x™ = h} (s), s = 8°,). Calculating the derivative with 
respect to x at point A? where Fo(x(s°,),u,s°,) = 0, we 
have 


ð Fo (yus) 
T 


=n (sj g =s. 





_ OGg(z,u,s) 


Fe Ba (hm(5), 8m) — 


m 





=h? (s), s=59, 
-h2 (s2,) on (GS. Sa) = Qom U): 
Calculating the derivative with respect to s at the same 
point, and taking into account the identity 
OP, 
ðs 





= r9o(x, 8), 


we obtain 


Ə Fo (£x,u,s) 
ðs 


z=h9, (5), s=89, 


= Gg (hin (Sm)> Us 8m) Be (hin ($8), 8m) = Pom (wu): 





All functions that enter the expression for OFo(2,t8) are 


continuous in the vicinity of Ap and P,,(u) is a finite 
quantity; therefore, OF o(e,th8) is continuous at Aj’. Thus, 
if a given u = koav/e — 1 satisfies the conditions 


u = koav'e — 1 > h? (s), 


where, for a given s > 1, h? (s) is the mth root of the 
function 


Qom (u) £ 0, 


P(x, s) = Jo(sx)Y, (x) — Jı (£)Yo (sx), 


the Cauchy problem (32) is uniquely solvable. 

Thus, we have demonstrated that, for arbitrary s = b/a > 
1 that may be sufficiently close to unity, there exists a root 
xt = x! (s) of DE (16) [more specifically, implicit function 
xt (s) specified by equation (16) in the vicinity of point AQ]. 
In other words, for an arbitrarily thin dielectric layer there 
exists (at least one) surface wave propagating in the Goubau 
line under the condition 


u = koav'e — 1 > h(s). (34) 


Relationship (34) is the necessary condition that provides the 
existence of at least one such wave. Taking into account the 
domain of function Fo(x,u,s) and Fo(x,u,s) we have the 
condition x € (0, koa,/e). 


V. CROSS-SECTIONAL STRUCTURE OF HIGHER-ORDER 
WAVES 


The cross-sectional structure of higher-order waves, or 
numbers of zeros (oscillations) of potential function ¢(7), can 
be determined by using the classical Sturm—Liouville oscilla- 
tion theory which may be a subject of a future separate study. 
In this paper, we prove the necessary mathematical items 
concerning the structure of higher-order waves independently 
using the analysis based on explicit expressions of the 
functions entering the DE. The result is: if {x,}/£, € [0, u] 
are K real zeros of function (16) w.r.t. x, then K = n 
where n = 1,2... is the sequence index of potential function 
én(r). For À going to infinity, the number of zeros of a on (r) 
defined on (a,b) also goes to infinity, while the intervals 
between two neighboring zeros contained in the interval tend 
to zero. 

Figure 7 illustrates this behavior for the first four potential 
functions. We see that the potential function corresponding 
to the principal wave (superscript zero) has no oscillations on 
the interval covering the cross section of the dielectric layer, 
while the potential functions corresponding to higher-order 
waves (superscripts 1,2, and 3) have, respectively, 1,2, and 3 
oscillations on this interval. 
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Getx,u,1.73)10 
Fa(x, 1.78) 


0.3 
0.2 
P1(1.78 ,x) 


0.1 


0.0 





-0.1 











Fig. 5: Graphs show the location of the first roots of DE (16) w.r.t x for 
u = 4 and s = s9 = 1.78: the first three roots are close to the zeros of 
function ®2 corresponding to vertical asymptotes with u = hess) where 
Py (Sy. u) = 0. The root ro (s9) of DE (the first zero of Fo) corresponds 
to the principal surface wave and is close to the first zero of ®2(z, s9) 
(the Ist vertical asymptote of function Fy), the second root xı(s9) = u 
corresponds to the first higher-order wave and coincides with the first zero 
h? (89) of function ®1. 
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Fig. 6: Location of the first four roots of DE (16) w.r.t x for u = 4 and 
s = 3.4: the roots are close to the zeros of function ®2, corresponding to 
vertical asymptotes on the graph. 


VI. RESULTS OF CALCULATION OF SURFACE WAVES IN 
THE GOUBAU LINE 


By solving numerically the Cauchy problems (28) and (32) 
we calculate the normalized propagation constants y of the 
fundamental (with no oscillations) and higher-order surface 
modes of the Goubau line vs. parameter s for different values 
of parameter u. The number of zeros of potential functions 
d(r) corresponds to the index of a higher-order wave; the 
cutoff values of s are determined by condition (2). The results 
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Fig. 7: Potential function ¢(r) (10) of the fundamental and higher-order 
surface waves of the Goubau line calculated for different eigenvalue indices 
n at a = 2, b = 6.8, and € = 5. 


are summarized in Fig. 8; it can be seen that the potential 
(eigen) functions preserve their types (number of zeros or 
oscillations) along the curves on the phase plane (s, y). 

Tables 1 and 2 exemplify some numerical values of 
the normalized propagation constants of the GL surface 
modes calculated for normalized (dimensionless) values of 
the geometrical parameters. The number of decimals show 
the accuracy of calculation using numerical solution of the 
Cauchy problems. 


Tab. 1: Roots x of DE (29) and the corresponding values of y at a = koa = 
0.2 and € = 2 


ay b] e] s] x] 7 
02 T.0060 
02 1.0133 
02 1.0201 
02 1.0277 
02 1.0347 
02 1.0417 
02 1.0290 
02 1.0562 
02 1.0632 
02 1.0702 


Tab. 2: Roots x of DE (29) and the corresponding values of y ata = 1 
and € = 2 


T oe a el 
j ity 2f [_09882 [10107 
| 12) 2f 12) oom [_1.0280 
| 13, 2f 13) 09506 | 10m 
| 14) 2f 14| 0926| 104 
| is) 2| 15| 03% | 107 
| 16) 2] 16| 0869| LIS 
E i7 2f 17| osso 1.1392 
t| 18| 2f 18| 08046 | 11630 
I| 19| 2f 19| o7 [1.1860 
1| 20| 2| 20) 0762| 1275 
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Fig. 8: Normalized propagation constants ~y (vertical axis) of the fundamental (with no oscillations) and higher-order surface modes of the Goubau line v.s. 
s calculated for u = u1 = 8 (blue), ug = 4 (black), and u3 = 2 (red). Inserts show potential functions (r); the number of zeros corresponds to the index 


of a higher-order wave. 


VII. CONCLUSION 


In this work, we have shown that a principal surface mode 
in the Goubau line exists for arbitrary values of permittivity 
of the cladding and radius of the inferred conductor (wire). 

If the electrical cross-sectional dimension of the line 
is sufficiently large, koave— 1 > hY(1) where h{(1) is 
an absolute constant, then there exist several higher-order 
GL surface modes, each having a distinct cross-sectional 
structure in the dielectric layer; namely, the number of zeros 
(oscillations) of the potential function corresponding to the 
principal and the nth higher-order GL mode equals the index 
of the wave: n for the principal and n for a higher-order 
mode. 

The longitudinal wavenumbers of the higher-order modes 
can be efficiently calculated by deriving the numerical so- 
lution of a Cauchy problem obtained by the parameter 
differentiation method. 

We have also reduced the determination of the spectrum 
of symmetric surface modes to singular Sturm-Liouville 
boundary eigenvalue problems on the half-line. The latter is 
reduced to the solution of (16) DEs involving the cylindrical 
tangent and cotangent functions, and then to the determi- 
nation of the parameter dependence of the DE zeros using 
analytical and numerical solution to the Cauchy problems. 

Several limiting cases have been investigated. We show 


that when an infinitely thin, perfectly conducting rod is 
placed into the cross-sectional origin of a dielectric fiber, 
it creates a regular perturbation of the spectrum of surface 
waves of the Goubau line with respect to the surface wave 
spectrum of the fiber. The spectrum of surface waves of the 
Goubau line with an infinitely thin dielectric coating is a 
perturbed set of zeros of a well-defined family of functions. 

The present approach remains valid for a GL with a lossy 
dielectric layer (complex permittivity parameter). In that case 
it is necessary to investigate all functions entering the DEs 
as functions of one or two complex variables and apply the 
methods of searching for complex roots of the DEs. This is 
a subject of our future study. 
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